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, | ■ We investigate excitation of Kaluza-Klein modes due to the parametric resonance caused by 

O |' oscillation of radius of compactification. We consider a gravitational perturbation around a D- 

, dimensional spacetime, which we compactify on a (D — 4)-sphere to obtain a 4-dimensional theory. 

The perturbation includes the so-called Kaluza-Klein modes, which are massive in 4-dimension, as 
well as zero modes, which is massless in 4-dimension. These modes appear as scalar, vector and 
second-rank symmetric tensor fields in the 4-dimensional theory. Since Kaluza-Klein modes are 
troublesome in cosmology, quanta of these Kaluza-Klein modes should not be excited abundantly. 
However, if radius of compactification oscillates, then masses of Kaluza-Klein modes also oscillate 
qq . and, thus, parametric resonance of Kaluza-Klein modes may occur to excite their quanta. In this pa- 

per we consider part of Kaluza-Klein modes which correspond to massive scalar fields in 4-dimension 
and investigate whether quanta of these modes are excited or not in the so called narrow resonance 
£SJ ■ regime of the parametric resonance. We conclude that at least in the narrow resonance regime 

quanta of these modes are not excited so catastrophically. 

On " 
On 

"o 
cr 



I. INTRODUCTION 



. £^ \ When we consider unified theories, spacetime may have more than four dimension. For example, the critical 
dimension of the superstring theories [jl| is equal to ten, and the low-energy effective theory of M-theory Q is the 
eleven-dimensional supergravity. On the other hand, we recognize only four dimensional spacetime. Therefore, in 
order to describe our universe we have to adopt some mechanism of dimensional reduction. Compactification of extra 
dimensions on some compact manifold (the Kaluza-Klein prescription |||) has been conventionally adopted, while an 
alternative prescription was recently proposed by Randall and Sundrum [Q. Since the Randall-Sundrum prescription 
cannot be applied directly to more than 5-dimcnsion, in this paper we consider the Kaluza-Klein prescription. 

However, there are too many ways of compactification to specify a 4-dimensional theory from the higher dimensional 
unified theories. As for the superstring theory, although it is usually compactified on a Calabi-Yau manifold, there 
are so many Calabi-Yau manifolds. 

Fortunately, in some cases it may be possible to judge which way of compactification is acceptable since different 
compactifications may give different 4-dimensional theories and, thus, different cosmologies. In particular, properties, 
eg. mass spectrum, of the so-called Kaluza-Klein modes reflect the manifold on which spacetime is compactified. 
Thus, it may be effective to investigate effects of Kaluza-Klein modes on the 4-dimensional cosmology. 

In this respect, Kolb & Slansky Q considered the effects of Kaluza-Klein modes on our 4-dimensional universe. For 
simplicity, they considered a scalar field in 5-dimensional spacetime of the form M 4 x S 1 , where M 4 is a 4-dimensional 
FRW universe and S* 1 is a circle. They investigated Kaluza-Klein modes of the scalar field and showed that the Kaluza- 
Klein modes are dangerous in cosmology. To be precise, they showed that, if energy density of a Kaluza-Klein mode 
is not negligible compared with radiation density at an early epoch, then the Kaluza-Klein mode eventually dominate 
the universe. Thus, quanta of Kaluza-Klein modes should not be excited abundantly. It is expected that this result 
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holds in quite general situation since their arguments are based on the momentum conservation along the compact 
direction. Thus, quanta of Kaluza-Klein modes of any fields should not be excited abundantly. 

However, there is a possibility that quanta of Kaluza-Klein modes may be excited. First, the radius of compactifi- 
cation appears as a scalar field in 4-dimension called a radion Hence, if the radion potential has a local minimum 
then it can oscillate around the local minimum. The oscillation of the radion leads to oscillation of masses of Kaluza- 
Klein modes since the masses are proportional to the inverse of the compactification radius. On the other hand, 
it is well-known that oscillation of mass leads to the so-called parametric resonance phenomenon |^||. Hence, it is 
expected that, if the radion potential has a local minimum then quanta of Kaluza-Klein modes may be excited by 
parametric resonance JTo| |. 

Therefore, it is possible to obtain some knowledge about a way of compactification by investigating whether quanta 
of Kaluza-Klein modes are excited by the parametric resonance due to the radion oscillation. If the quanta are 



produced abundantly then the way of compactification should be rejected. In Rcf. |10|, the parametric resonance of 
Kaluza-Klein modes of a scalar field was investigated in the regime where the resonance band is narrow (the narrow 
resonance regime). It was shown that, for the case of compactification by a d-dimensional sphere, quanta of the 
Kaluza-Klein modes are not overproduced. 

However, since there should be many other fields (eg. gravitational fields, anti-symmetric fields, spinor fields, and 
so on), the investigation of only a scalar field is not enough: in principle, we have to investigate Kaluza-Klein modes 
of all fields in the theory. In this paper, we investigate a gravitational field. 

This paper is organized as follows. In section [Il| we describe our model and derive the action for our system. 
In section III we investigate whether catastrophic creation of quanta of the Kaluza-Klein modes occur or not. In 
section 



IV we summarize this paper. 



II. MODEL CONSTRUCTION 

In this s ection we derive an effective action for Kaluza-Klein gravitational modes as well as for background fields. In 



subsection II A we describe our model and basic assumptions verbally. In subsection [IB we calculate the tree action 



of our system from the Einstcin-Hilbert action. In subsection II C we stabilize the radion potential by the so-called 
Casimir effect to obtain the effective action for our system. 

A. Basic assumptions 

As already explained in Sec. ||, our purpose in this paper is to judge whether a particular model of dimensional 
reduction is acceptable or not by investigating parametric excitation of Kaluza-Klein modes. For this purpose, we 
have to specify higher dimensional theory in .D-dimension somewhat. Moreover, in order to obtain the 4-dimensional 
effective theory from the D-dimensional theory, we have to specify a particular model of dimensional reduction. 
Namely, we consider the following setting. 

(a) Let us consider a D-dimensional theory which includes gravity, N s scalar fields, Nd Dirac fields, and other fields. 
The gravity, scalar fields and spinor fields are assumed to be described by the Einstein-Hilbert action with a 
cosmological constant, the Klein-Gordon action and the Dirac action, respectively. There may be interaction 
among these fields. 

(b) We adopt the conventional Kaluza-Klein description of compactification: we compactify the D-dimensional 
spacetime on a d-dimensional compact manifold (d — D — 4) to obtain the 4-dimensional spacetime. As the 
compact manifold we take a (i-dimensional sphere, which we shall denote by S . 

(c) We consider the so-called Casimir effect to stabilize a radion potential, where the radion is a 4-dimensional scalar 
field corresponding to c ompa ctification radius In order to stabilize the radion potential by the Casimir effect, 



as shown in subsection II C, we have to assume that d > 2. 



(d) We also consider perturbations of D-dimensional fields around the background given by the Kaluza-Klein pre- 
scription. These perturbations are described as 4-dimensional massive fields called Kaluza-Klein modes. 

For this situation, in order to make calculations possible, we assume the following assumptions. (See the beginning 



of section III for further assumptions.) 



(i) The number N s (or Nd) of scalar fields (or Dirac fields) is sufficiently large compared with the number of other 
fields included in the ZJ-dimensional theory. 
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(ii) The cosmological constant induced in 4-dimension is zero. 

(iii) Terms in the 4-dimensional effective action which are order of 0(ip 3 ) can be neglected, where ip denotes Kaluza- 
Klein modes. 

The assumption (i) makes it possible to calculate corrections to the radion potential due to the Casimir effect. In 
principle, all fields included in the Z3-dimensional theory should contribute to the Casimir effect. However, if N s or 
Nd is large enough then contributions from other fields are small compared with those from the scalar fields and the 
Dirac fields. Hence, in this limit, the correction to the radion potential is parameterized by only N s and Nd- In fact, 
we can show that, after redefining the ZJ-dimensional cosmological constant, the corrected potential is parameterized 
by only two constants: the _D-dimensional cosmological constant and a particular combination of N s and Nd- 

Therefore, the assumption (i) combined with the assumption (ii) det ermines the corrected radion potential uniquely 
up to only one constant. Hence, as we shall show in subsection |H Cj explicitly, the potential is parameterized only 
by the present value bo of the compactification radius. Since bo can be eliminated from all equations of Kaluza-Klein 
modes by redefinition of variables, all information about the corrected radion potential can be obtained. 

On the other hand, the assumption (iii) makes it possible to investigate each Kaluza-Klein mode independently, 
provided that kinetic terms and mass terms are properly diagonalized. In this paper, we consider only a particular 
class of Kaluza-Klein modes corresponding to D-dimensional gravitational perturbations. (In Ref. [|lCfl a class of 
Kaluza-Klein modes corresponding to a D-dimensional scalar field was considered.) It is the assumption (iii) that 
makes it possible to consider these modes independently. 



B. Perturbed Einstein-Hilbert action 



We consider the D-dimensional Einstein-Hilbert action with a cosmological constant: 



I EH - 

2n z 



d u x*J=§{R-2k), (1) 



where n is a positive constant, R is the D-dimensional Ricci scalar, and A is the cosmological constant. We consider 
a gravitational perturbation Hmn around a background metric 9m ni which we shall specify below. 

9mn — s'mn + h-MN- (2) 
By substituting Eq.(||) into Eq.(Q), we obtain the perturbed Einstein-Hilbert action as follows. 

Ibb = ± j d°x^¥» [i? (0) - 2A - h™ (R$ n - + 



I (h* - 2h MN h MN ) m + \ (2h MM 'h N M , - hh MN ) R<g N 



\ {h M %, (2h% N h M f) + h., M {K> M - 2h M %) } 
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-a (hi z - \ii J ) fonn 



(3) 



where ";" denotes the covariant derivative compatible with g$ N , and R^n an< ^ are the Ricci tensor and scalar 
constructed from g^N- (^ ee Appendix |A| for a detailed derivation.) 

As for the background geometry, we compactify it on a d-dimensional sphere S d (d = D — 4) : 

g ( MNdx M dx N = g^dx^dx" + b 2 n ( i fdx l dx j , (4) 

where g^ and b are a 4-dimensional metric and a 4-dimensional scalar that depend only on the four-dimensional 
coordinate x^ (// = 0,1,2,3), and flij is a metric of the unit d-sphere which depends only on the coordinates x l 
{i = 4, • • • , D — 1). Hereafter, we denote Z?-dimensional, 4-dimensional and d-dimensional indexes by capital Latin 
(M, N, ■ ■ •), Greek (/i, v, ■ ■ ■) and small Latin ■ ■ •) letters, respectively. Note that the scalar b can be interpreted 
as radius of S d . 

Since it is convenient to analyze in the so called Einstein frame, we perform the following conformal transformation 
to make the new frame to be the Einstein frame. 
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9nu 



(5) 



Note that, if we take b to be the present value of b, then the conformal factor becomes unity after b settles to bo- 
Hence, the change of the frame does not affect the final result, provided that all physical results are interpreted after 
b settles to bo. On the other hand, our analysis becomes much easier if we adopt this conformal transformation. After 
the conformal transformation, the D-dimensional line element is written as 



9mn^ x dx 



(0) 



For this particular background, the gravitational perturbation Hmn can be expanded by harmonics on S d as follows 
h MN dx M dx N = [h l ™Y lm dx»dx» + 2{h\^(V (T)lm ) l + h$ )M (V {L)lm ) i }dxf>dx i 

(7) 

where Yj m is the scalar harmonics function; V(T); m & n d V(£v m are the vector harmonics; T(T)Zmj T(LT)imi an d TrLL)lm 
are the tensor harmonics. Here, the coefficients /ijjJJ, fc, ^(iW Mt)' ^(™t)> ^(ll) an< ^ ^(™) depend only on the 
four-dimensional coordinates x M , while the harmonics depend only on the coordinates x % on S d . (See Appendix [b| for 
definitions and properties of these harmonics.) 

Although this expression of Hmn includes many terms, some of them represent degrees of freedom of coordinate 
transformations. In fact, it is shown in Appendix O that, after gauge-fixing and redefining g^ u and b, the perturbation 
fiMN can be expressed as follows. 

)Y lm dx»dx v + 2h i ^{V (T)lm ) l dx»dx l + {h l ffi(T( T ) lm )ij + h l fo(T (Y)lm ) ij }dx i dx> 



h M Ndx M dx N = ^ \h l ™Y •'• -''•'<" •>'■'"' n ~ 1 1 ' T 1 



(8) 



where the summations are taken over 2 > 1 for the scalar and vector harmonics, and over I > 2 for the tensor 
harmonics. 

Finally, by substituting this expression into Eq. (^), we obtain the following action. 



Ieh = I (0) + + I {2) + 0(h% 



where 



j(o) 

j(2) 



2 k 



E/^v^)(4r+4r+4r), 

Lm 



(9) 

(10) 
(11) 



and is linear in ft,. When a total action of the system is considered (see the setting (a) in subsection II A), 1^ 
should be canceled by other linear terms in the total action because of the equations of motion. Hence, there is no 
need to consider The constant k, the scalar field a, the potential Uo, and Lagrangian densities C 



(T,V,Y) 



defined as follows. First, the constant k (> 0) is defined by 



Next, a is a scalar field defined by 



a = ( r ln(- 



Co 



d(d + 2) 



2k 2 



are 
(12) 

(13) 
(14) 
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and is called a radion field M. The potential Uo(cr) of a is given by 



AC 



d(d-l) -(d+2)^- 
K-e ff o . 

2k 2 6 



(15) 



Finally, the Lagrangian densities £j 



,(T) 



(T,V,F) 



are given by 



where 



/•P0 _ /»0^)r fjm rim 
J -lm — Hm I ' h tkv ) n (Y) 



(X) lm lm 
X X j 



(16) 
(17) 
(18) 



A 



(T) 



+e 



lm 



b - 



2 d+2 7 



.2 "(T) 



l) + (d- 


l)(d- 


2)}&o" 2 


(3- 


jd(d + 










(/ > 2). 







2e~ d ^A 



(19) 
(20) 



(Since hereafter we analyze M™% (or x irra 



only, we have not written down explicit form of and £jn? - ) Hence, 



up to the second order in Kaluza-Klein modes, h$\ (or x' m 



is decoupled from all other Kaluza-Klein modes. (See 

the assumption (iii) in subsection II A.) In this paper, for simplicity, we investigate th e par ametric resonance of x lm 
(or /i|™-)) only. It should be stressed again that it is the assumption (iii) in subsection II A that makes it possible to 
consider these modes independently. 

Here, note that the constant &o is arbitrary at the moment and that we can eliminate it from all equations. However, 
in the next subsection, we chose bo to be the present value of b. 



C. Stabilization of the radion potential 

Here we assume that the Casimir effect [Q stabilizes the radion potential. In other words, we ass ume that quantum 
corrections to Uq stabilize compactification. Thanks to the assumptions (i) and (ii) in subsection II A , we can easily 
calculate the stabilized radion potential, which we shall denote by U\. 

First, for all Kaluza-Klein modes of N s scalar fields and iVj Dirac fields, mass in the Einstein frame is proportional 
to the common quantity. 

M 2 oc bo 2 e-^ d+2 ^^°. (21) 

Hence, because of the assumption (i), we can obtain the following 1-loop effective potential. (See, for example, Ref. |1(J 
for details.) 

V lloop (a)=Ce-^ d+2 ^, (22) 

where C is a constant determined by N s and N^- There are corrections to the Z?-dimensional cosmological constant 
A, too. However, these corrections can be absorbed in A by redefining A. Hence, the corrected radion potential U\ is 
given by 

loop {p )• (23) 
Next, the assumption (ii) requires that the minimum of U\ should be zero. Hence, A is determined as follows. 

A - d(d-l)(d + 2) 

A ~ 2(d + 4)b 2 ' (24j 

where we have specified the constant fro as 
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' d(d-l) 
(d + 4)Ck 2 



Thus, we obtain the following radion potential Jll| , [l0| . 

Ui(a) 2 
where a is a constant given by 



-2(d+2)a/a 



-da /(To 



d + 2 



d+2 



-{d+2)a/a 



a 



d(d- l)(d + 2) 
2(d + 4)n 2 b Q 2 ' 



(25) 



(26) 



(27) 



Thus, in this case, the 4-dimensional 
1 then the potential is completely flat. 



This potential has the minimum zero at a — 0, provided that d > 2 
cosmological constant is actually zero and &o is the present value of b. If d 
Hence, hereafter, we assume that d > 2. 

It is notable that the radion potential has been determined uniquely up to the present value bg of the compactification 
radius b. Since bo can be eliminated from all equations of x lm by redefinition of variables, all information about the 
corrected radion potential has been obtained. 

Finally, the perturbed Einstein-Hilbert action corrected by the Casimir effect is given by Eq. (|J) , provided that Uq 
is replaced by U\ given by Eq. (26). Therefore, hereafter, we investigate the dynamics of the field x lm coupled to the 
radion field a in the expanding universe. These fields are described by the following action. 



(28) 



where 



d*xy/-g(°) 



(29) 



Here, we mention that in the total action there may be other terms which are second order in \ lm . In fact, there 
should appear other mass terms for \ lm since some zero modes may couple to x' m - However, providing the assumption 
(vii) in the next section, these additional mass terms are small enough since zero modes are considered as usual low- 
energy matter whose energy density and pressure are bounded from above by 4-dimensional cosmological parameters. 
Therefore, the above action is enough for our purpose. 



III. PARAMETRIC RESONANCE 

In this section we investigate parametric resonance of the field x lm due to the oscillation of the radion field a in 
the expanding universe. For this purpose, we specify 4-dimensional geometry and initial condition somewhat. 

(iv) The background 4-dimensional geometry is expressed by the flat FRW metric, and the radion field is homoge- 
neous in the FRW universe. 

(v) Initially, there is no excitation of Kaluza-Klein modes. Thus, the FRW universe is driven by the radion field as 
well as by usual homogeneous matter corresponding to zero modes. 

Since the radion potential Ui(a) has exponential terms, it seems difficult to analyze oscillation of a with large 
amplitude analytically. Hence, we assume that the amplitude of the oscillation is small enough. Moreover, for further 
simplicity, we assume that the time scale of the cosmological expansion is much longer than the time scale of the 
radion oscillation. Namely, we assume the following. 

(vi) We assume that the amplitude of the radion oscillation around a = is small compared with <jq • 

(vii) The energy scale of the cosmological expansion is much lower than that of compactification: we assume that 
Hbo <C 1 and Hb^ <C 1, where H = a /a and the dot denotes the derivative with respect to the cosmological 
time. 

Note that the assumption (vi) leads to the so called narrow resonance regime for the parametric resonance. 
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A. FRW background and homogeneous radion 



According to the assumptions (iv) and (v), we consider the homogeneous radion field and the flat FRW universe, 
which is driven by the radion field and usual matter corresponding to zero modes: 



g^dx^dx" = -dt 2 + a(t) (dx 2 + dy 2 + dz 2 ), 



a = a(t). 



(30) 



The cosmological evolution equations are 

H 2 = 
H + H 2 



8irG 



N 



3 

4nG N 



(Pa- + Po) 

[{Pa + Po) +3(p a +Po)} 



(31) 



where the subscripts 'c' and '0' mean contributions from the radion a and those from usual matter corresponding to 
zero modes, respectively. Here, H = a/a, Gn — k 2 /8tt, and the dot denotes derivative with respect to the time t. 

Now, because of the symmetry of the background, it is convenient to analyze the field \ lm in the momentum space. 
Hence, we decompose the field \ lm as follows: 



Xim = a 



(32) 



(12) 

where x ; 'u are re& l functions of the cosmological time t and the vector x represents (x,y,z). Correspondingly, the 



action I x is rewritten as 



J v = / dtL 



iv = 



l,m i=l,2 



where 



a 2 I 9 2 \a 



(33) 



(34) 



Although the kinetic term in the above Lagrangian has time dependence, we can eliminate it by using a new time 
variable r defined by 



dr 

dt =GXP 



Co 



(1 2) 

In fact, the equations of motion for are written as 

d 2n 



Or 2 



+ n 2 (T)Q = 0, 



(12) 

where Q denotes x ; ' ^ and Q is defined by 

Q 2 ( T ) = m 2 (t(r))exp 
The hamiltonian with respect to the time r is 

H = t(P 2 + n 2 Q 2 ) 

where P is the momentum conjugate to Q. 

After the oscillation of a ends up at a — 0, the right hand side of Eq. ([55|) becomes unity. Thus, after the end 
of the oscillation, the hamiltonian Ti with respect to r coincides with the hamiltonian with respect to the original 
cosmological time t. At the same time, those hamiltonians become conserved quantities. Hence, we can analyze the 

(1 2) 

excitation of % ; ' ^ by using the new time variable r. 



(TO 



(35) 



(36) 



(37) 



(38) 
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B. Small oscillation of the radion 



Providing the assumptions (vi) and (vii) , the equation of motion for a gives 



^ =5-(f)cosw(t-f ), (39) 



where to is a constant, to is defined by 



"=\fim=y d ^ (40) 



1 ^fi-^^-in). 

_ + 4) {1(1 + d- 1) -T ■ <<• - - I!.) - Uj - - J I!./ -:- ~ , 

6g ^2 . ^ 



and (<gC 1) is a slowly varying amplitude: 

er(t) oc a~ 3/2 . (41) 
Now, by using Eq. (|3^), we can express f£ 2 in terms of the new time variable r up to the first order in a: 

n 2 {r) = lu 2 [A + B + ecos{2w(r - r )} + 0{e 2 )] , (42) 

where tq is a constant and 

2tl(^M 2 /9 2 3 - 

2(d + 4) {?(/ + d - 1) + (d - 2){d - 1)} - d(d - l)(d + 2) 

e = aUd - l)^ + 2)(d + 8) -^^ 1)( ^ 2)(d+10) - A (fcVa)2 + 4(d + 4)1 . (45) 

The first term in A is essentially the squared momentum along the three-dimensional spatial directions divided by 
uj 2 , B is essentially the squared momentum along S d divided by u> 2 , and e is the properly normalized amplitude of 
the radion oscillation. 

C. Result 

The expression ([l2]) is of almost the same form as the expression of tt 2 investigated in Ref. poft . Hence, we can 
analyze the parametric resonance in the same way as in Ref. juj. Thus, the necessary and sufficient condition for the 
catastrophic creation of quanta of the Kaluza-Klein mode is 

B<1, 

\B-l\=0(e). (46) 

However, it is easily shown that 

B > 3 (47) 



for any values of d (> 2) and I (> 2). (See the setting (c) in subsection HA.) Therefore, it is concluded that quanta 
of the field x lm are not overproduced by the parametric resonance due to the small oscillation of the radion a. 

IV. SUMMARY AND DISCUSSIONS 

We have investigated excitation of Kaluza-Klein modes due to the parametric resonance caused by oscillation 
of radius of compactification in the narrow resonance regime. In particular, we have considered a gravitational 
perturbation around a Z?-dimensional spacetime compactified on a (D — 4)-sphere. Among many Kaluza-Klein modes 
of the gravitational perturbation we have investigated part of Kaluza-Klein modes which correspond to massive scalar 
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fields in 4-dimension. We have obtained that quanta of these modes are not excited so catastrophically in the narrow 
resonance regime. 

Now, several comments are in order. 

First, we discuss about the consistency among the assumption (v), (vi) and (vii). Energy density p G and pressure 
p a of the radion is of the following order. 

Pa -Pa- €% 2 K-\ (48) 

where e ~ <t/<7q. Thus, the condition H 2 ~ n 2 p a can be satisfied, even if e <C 1, Hbo <C 1 and Hb^ <C 1. In other 
words, the three assumptions (v), (vi) and (vii) can be satisfied simultaneously. 

Next, we discuss about validity of our strategy. Although we have investigated only a particular class of Kaluza- 
Klcin modes in this paper, in principle we have to investigate all Kaluza-Klein modes in order to judge whether 
a way of compactification is acceptable or not. If we find that one of these many Kaluza-Klein modes is excited 
abundantly, then the way of compactification should be rejected. In this respect, there may be an embarrassing 
expectation: it is easily expected that, even if we find that Nkk 3> 1, the result will be significantly changed by 
backreaction, where Nkk is the number of created quanta of a Kaluza-Klein mode. In fact, by the energy conservation 
and the assumptions (v) and (vii), it is easily concluded that Nkk *C 1. Thus, it might be expected that under the 
assumptions (v) and (vii) we might not be able to judge whether a way of compactification is acceptable or not. 
However, this embarrassing consideration does not spoil our analysis. Actually, if we obtain Nkk S> 1 by ignoring 
the backreaction, then this overestimating result implies at least that pkk is not negligible compared with p a since 
the backreaction becomes effective only when Pkk is not negligible compared with p a , where pkk is energy density 
of the Kaluza-Klein mode. Therefore, providing the assumption (v), the result Nkk S> 1 obtained by ignoring the 
backreaction implies that pkk is not negligible compared with energy density of radiation, which should be rejected 
as explained in Sec. [|. Because of almost the same reason, we can justify the assumption (iii). First, provided that 
the higher-order coupling constant is at most of order unity in the unit of b$, those higher-order terms should be 
inefficient as long as Pkk <C &^ 4 . Physically, this statement is almost equivalent to the statement that sufficiently 
small oscillation does not feel any higher-order corrections to its potential. Therefore, we can neglect higher-order 
terms if the created number of Kaluza-Klein quanta is small. (See the assumption (v) for the initial condition.) In 
other words, if higher-order terms becomes efficient at some epoch, then it implies that the number of Kaluza-Klein 
quanta created so far is not small, and such a situation should be rejected as explained in Sec. |l]. 

Finally, we comment on possible future works. (A) We have to investigate Kaluza-Klein modes of other fields as 
well as other gravitational modes. (B) We have to take inflation into consideration. (C) If we can extend our analysis 
to the broad resonance regime, then it will give stronger constraints on ways of compactification. Note that, even 
in the broad resonance regime, Eqs. (|36| ) and ( |38| ) can be used. (D) Other mechanisms of compactification might 
deserve considering, for example, the monopole-like configuration of an anti-symmetric field [T^ ]. (E) It might be also 
interesting to investigate the parametric resonance of Kaluza-Klein modes in the Randall-Sundrum scenario [EJ with 
stabilized radion |l3| . 
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APPENDIX A: SECOND-ORDER EINSTEIN ACTION 

The purpose of this appendix is to expand the Einstein-Hilbert action by Hmn up to the second order. Although 
in the main part of this paper all quantities defined in /^-dimension have "~" , we omit it in this for simplicity. 
Let us decompose the metric gMN into the background g^j N and perturbation Hmn- 

9ain — 9mn + h-MN- (Al) 
Hereafter, we use the background metric g^N t° ra ise the indices of Iimn'- 
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h M . M _ „MOl 

/i MJV = fl (°) M V 0)JVP /iop, 

fc^^W (A2) 

First, we can easily expand g MN , v /= 5 an d the Christoffel symbol rj^-jy as follows, where g is the determinant of 

9M N- 

g MN = g MN (0) _ h MN + h MO hQ N + 

^—g = {l + \h + \ (\h 2 - \ h MN h MN ) } + 0{h% 

^mn = r {a) MN + r (1) M^ + r {2)P MN + o(h 3 ), (A3) 

where and T^ MN arc tnc determinant and the Christoffel symbol of the background metric g M \, and 

F^MiV = 2 ^ {h-QM , N + ft-QW , M - hMN , q), 

MAT = — 2 h PQ {flQM , N + flQN , M — h-MN , q)- (A4) 

Here, " , " denotes the partial derivative. 

Next, the Ricci tensor is expanded as follows. 

Rmn — Rmpn — r£f N i P — r£f P _ ^ + t p p t mn — t p n t mp 

= #MiV + R MN + ^MiV + °(^ 3 ): (A5) 

where JJ^ is the Ricci tensor for the background metric g M ^ Nl and 

^mjv = T^mn , p — T^mp , at i (A6) 

p(2) _ r (2) P r (2) P , r (l) P r (l) S r (l) P r (l) S / a y\ 

H MN - 1 MN , P - 1 MP , AT + 1 SP 1 MAT ~ 1 SAT 1 MP- ( A <) 

Therefore, we obtain the following expansion. 

V~9R = V~99 MN Rmn = (^ (0) + R {1) + R (2) + 0(h 3 )) , (A8) 

where R^ is the Ricci scalar for the background metric g M ^ N , and 



= f - W 0) 9M N ) + (h MN ; N - ^ M );M, 



2 



i? (2) = ^ - v^pq) i? (o) + ^ft/e + 9 w mn rw 

_I ^MATpCO) _.MJVn(l) i ^JPrW 
2 n MN n n MN ' 2 " MN 

MATnd(O) 



MN 



l(h 2 - 2h p Qh PQ )R^ + ±(2/^^" - M MJV )# 



MN 



+ \ {h MN ;P + (2h P M , N ~ hMN' P ) + h, M M - 2/> MJV ; at)} + X P p . (A9) 

Here, X p is a vector constructed from hMN, and " ; " denotes the covariant derivative compatible with the background 
metric g$ N . 

Finally, by using the above expressions and integrating by part, we obtain the perturbed Einstein-Hilbert action 
up to the second order in hMN- 
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i EH = ^LJ d D xV=g(R (Q) - 2A) 

= h I dDx ^^ [ R(0) - 2A - hMN - \*f>A + aA 

+ \ (: 



2h MM ' h», - hh MN 



(o) 

MAT 



Ln M\N a MN 



:M 



(1> 



A I 



2h 



MN \ \ 



\\\h*-\h MN h M N 



0{h 3 ) 



(A10) 



APPENDIX B: HARMONICS ON D-SPHERE 

In this Appendix we give definitions and basic properties of scalar, vector and tensor harmonics on a unit en- 
sphere Jl4| ] . Throughout this Appendix we will use the notation that fly is the metric of the unit d-sphere and D is 
the covariant derivative compatible with f2y. 

1. scalar harmonics 

The scalar harmonics is supposed to satisfy the following relations. 

D 2 Y lm + l(l + d-l)Y lm = (i>0), (Bl) 

(B2) 



2. vector harmonics 



First, in general, a vector field V on the unit d-sphere can be decomposed as 

Vi = V {T )i + dif, 

where / is a function and Vr T ) is a transverse vector field: 

£> i (V (T ))i = 0. 



(B3) 



(B4) 



Thus, the vector field Vi can be expanded by using the scalar harmonics Yi m and transverse vector harmonics 

V(t) im as 



(B5) 



where cfefs and are constants, and the transverse vector harmonics V(t) im (I > 1) is supposed to satisfy the 



following relations. 



# 2 V ( r) im + W + d - 1) - 1} Vfa , m - 0, 

= 0, 

d d xVU n ij {V(T) lm)i (V( T ) ('m')j = fa' <*n 



^From Eq. (B5), it is convenient to define longitudinal vector harmonics VtL)im by 

(V( L )l m )i = diYi m (/>1). 



(B6) 



(B7) 
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It is easily shown that the longitudinal vector harmonics satisfies the following properties. 
D 2 V (L) lm + {1(1 + d-l)-(d- 1)}V (L) lm = 0, 

0'(Vfajm)< = -l(l + d-l)Y lm , 

D[i(V(L) lm)]} = 0, 
J d d xVU O ij '(V(i) lm)i(V( L )l'm')j =l{l + d-l)S W $mm>, 

J d d xVn n ij (v (T ) i m UV {L) Vm ,) 5 = o. (B8) 



3. Tensor harmonics 

First, in general, a symmetric second-rank tensor field TV,- can be decomposed as 

T ij = T (T) ij + A Vj + Dj V, + , (B9) 
where / is a function, V% is a vector field and Tit) ij is a transverse traceless symmetric tensor field: 

(T(T))i = 0, 

= 0. (BIO) 

Thus, the tensor field can be expanded by using the vector harmonics V(r) im ancl V{L) imi ancl transverse vector 
harmonics 7m i m as 

T ij = ^2 [ C \t)( T (T) lm)ij + C l {LT) { D i ( V (T) lm)j + A (V( T ) / m )j} 



lm 



+ C (™L) {A (V ( L) lm)j + Dj (V (L) lm )i} + C^Y^-j , (Bll) 

where c^, c|™ T ^, c|™ L ^ and c|™^ are constants, and the transverse tensor harmonics T(r) i m (I > 2) is supposed to 
satisfy the following relations. 

D 2 T (T) lm + {1(1 + d - 1) - 2}T (T) lm = 0, 

(T(T) lm)\ = 0, 
D Z (T(t) lm)ij = 0, 

J d d xVfi fl"' ^"(T( T ) j m )«(T (T) = <y M /*mm'. (B12) 

From Eq. (Bll), it is convenient to define tensor harmonics T(lt) inn ^(ll) imi anc > Try) im by 

(T(LT) lm)ij = A(V(T) i m )j + Dj(V(x) lm)i, (I > 2) 

2 

(T(LL)lm)ij = A(V(L)(m)j + A (V{L) lm)i ~ -j^ijD (V(L)lrn)k, (I > 2) 

(T(r) im)ij = ^i]Y(Y) lm (I > 0). (B13) 

It is easily shown that these tensor harmonics satisfy the following properties. 
D 2 T [LT) lm + {1(1 + d - 1) - (d + 2)}T (LT) lm = 0, 

D l (T (LT) lm )ij = -{1(1 + d - 1) - d}(V (T) lm )j, 

(T { LT)lm)i = 0, (B14) 

AT (Li) im + {/(J + d - 1) - 2d}T {LL) lm = 0, 

D l (T {LL) Im ) y = - 2(<i ~ + d - 1) - d}(V (i) Iro )j, 

(r (LL)tol ) 1 , = 0, (B15) 
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and 

D 2 T (Y) lm + l(l + d-l)T {Y) lm = 0, 

D % (T(y) lm)ij = \V{L) lm)ji 

{T (Y )im)i = dY lm . (B16) 

It is also easy to show the following formulas of integral as well as the orthogonality between different types of 
tensor harmonics. 

d d xVU ST" Q" (T (LT) / m )y(T (iT ) Vm >) V j' = 2{l(l + d- 1) - d}5 u > 5 mm >, 

d d xVU SV j '(T {LL) l m )i 3 {T {LL) Vm ,) Vjl = Al{l + d -^ d - l \ l { l + d - 1) - d}8 W 6mm', 

d d xVU ST 4 Q n (T(y) lm)ij(T(Y) l'm')i'j' = d 5u> 5 mm > . (B17) 

4. Integration formulas 

Finally, we list up some other formulas of integration. 

d d xVn D l (V (T) lm )j A(V (T) i'm') j = W + d - 1) - 1}6 W 5 mm> , 

d d xVn d\v {t) lm )i d z (v {t} Vm ,y = -(d - i)s w s mm ,, 

d d xVTl D\T [T) i m y k D t (T (T) Vm ,) ik = {1(1 + d - 1) - 2}5 W L.', 
(fWn i m y k D t {T {Y) Vm ,) jk = d{l(l + d - 1) - 2}«„, 

d d xVfi D k (T(T) IrnY 3 Di(T(T) l'm')jk = ~dSll' S m m' ) 

(?Wfi £ fc (T (y) Im ) y A(T(y) i'mOifc =l{l + d- l)Sw <W (B18) 

APPENDIX C: GAUGE FIXING 

In this appendix we show that there is a particular choice of gauge in which the gravitational perturbation hj^N is 
expanded as Eq. (||). 

First, the infinitesimal gauge transformation of Hmn due to a vector field £ M is expressed as 

h-MN ^ h MN — £m-,N — £n-,M, (CI) 

where £m = 9mn£ N i an< ^ denotes the covariant derivative compatible with the background metric g M \- 
Next, we expand hMN and £m by harmonics on S d , which are defined in Appendix [B| as 

h MN dx M dx N = J2 [ti™Y lm dx»dx» + 2{h[™ )ll (V (T)lm ) i + h$ ) p{V (L)lm ) i }dx> i dx i 

lm 

(C2) 

and 
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i M dx M = [C Y ^dx^ + {ti l (T)(V(T)lm)i + $L)(V W lr n )i}dx i 



(C3) 



Here, the coefficients h 



lm ulm ulm ulm 
fiui (T)jU' (L)m' (T)> 



■..In 



Jr. 



Ky)i £m > mti anc ^ depend only on the four-dimensional 



l (LT)' "(LL)' *>p ' S(T) 



coordinates x M , while the harmonics depend only on the c oord inates a:* on S . 

Hence, by substituting the above expansions into Eq. (CI), we obtain the infinitesimal gauge transformation of 
each coefficient. 




b\ d/2 
bo 



-d/2 



4(L) 



bo 



2Z(Z • 



1) 



(C4) 



where V is the covariant derivative compatible with the 4-dimensional metric g^ u . 

lm 
i 5 



Therefore, by choosing £^ m , and as 



>(T) 

tin i 
?(L) 

tZm 



f lm 
>(L) 



&-^ ( ™ T) (/ > 2), 
b" 1 /^ (i>2), 



(LL) 

&0 
— h lm 



(L) M 



&0 



d/2 



4(L) 



(i>l), 



bo 



d/2 



(C5) 



we can always make a gauge transformation such that 



ulm 

n (LT) ~ 


-> 


(i > 2), 


h lm 
(LL) 


-> 


(1 > 2), 


ilm 
\L)^i 


-> 


(/ > 1), 


ilm 

n (Y) 


-> 0. 





(C6) 

Moreover, h ® v and h.9jL can be eliminated by redefining the 4-dimensional metric g^v (or g^v) and 4-dimensional 
scalar field b (or a). It can be easily shown that the gauge transformation of ft,™ and ft9^L is equivalent to the 
4-dimensional gauge transformation. 

Therefore, we obtain the following expansion of Hmn'- 

h MN dx M dx N = £ [h l ™Y lm dx»dx" + 2h l ff )tl (V (T)lm ) i dx fl dx i + {h L ff } (T {T)lm ) tJ + h\™ } (T {Y)lm ) tJ }dx l dx3] , (C7) 

Irn 

where the summations are taken over I > 1 for the scalar and vector harmonics, and over I > 2 for the tensor 
harmonics. 

The remaining gauge freedom is given by and The first one generates the gauge transformation of h^\^ 



as 



(^(T))/i — * (^(t))ai ~b 



u-lflm 



(C8) 
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Another remaining gauge freedom £;| u generates the usual 4-dimensional gauge transformation. 
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